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The charge form factor of the pion is calculated in lattice QCD. The non-perturbatively improved 
Sheikholeslami-Wohlert action is used together with the 0(a) improved vector current. Other 
choices for the current are examined. The form factor is extracted for pion masses from 970 MeV 
down to 360 MeV and for momentum transfers Q 2 < 2 GeV 2 . The mean square charge radius is 
extracted, compared to previous determinations and its extrapolation to lower masses discussed. 
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I. INTRODUCTION 

QCD without doubt is the correct microscopic theory 
describing all strong interactions, a fact that has been es- 
tablished mainly by impressive agreement between the- 
ory and experiment in the perturbative sector. Com- 
paratively few results were obtained in non-perturbative 
QCD, which deals with physics on the scale of Aqcd or 
the size of a hadron. It is therefore an obvious challenge 
to derive the internal structure of a hadron from first 
principles, entirely within QCD Next to the nucleon, 
the pion is an obvious candidate for such an attempt. 
Its global features, like charge, spin and isospin, repre- 
sent no challenge and are trivially included in any model. 
Specific features that are actually testing details of our 
theoretical understanding are observables like the pion 
form factor or its polarizability. In this paper, we report 
on an extensive study of the pion form factor based on 
lattice QCD. First results were already reported in 0,O|. 

At first glance, the pion looks like a manageable two- 
body system and there have been many descriptions of 
the pion based on effective models or QCD inspired ap- 
proaches. One feature all these attempts share is that 
confinement, the most striking feature of QCD, is - in 
one way or another - put in by hand. This is of course 
an unwanted step when one sets out to calculate the pion 
form factor or its mean square charge radius, which re- 
flect the form and size of QCD confinement. Here one 
obviously wants to proceed from first principles, from 
the fundamental QCD Lagrangian itself. 
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Several papers have already dealt with aspects of the 
pion structure in lattice QCD. An often considered quan- 
tity is the 'Bethe-Salpeter amplitude', the relative quark- 
antiquark wave function, extracted from 2-point func- 
tions 0, 0, IE 0- Another approach has been based 
on gauge invariant density-density correlations Q, 0, 13 , 
most recently by Alexandrou et al. [9j . Two groups have 
already calculated the pion charge form factor, proceed- 
ing, as in the present work, in the quenched approxi- 
mation. The pioneering work was done by Martinelli et 
al. |lOl| . fo llowed by more detailed calculations of Draper 
et al. 11]. One of the findings of the latter work was 
that the pion form factor could be described quite well 
by a monopole form, as suggested by vector meson dom- 
inance. As shown in |2(, the range parameter is in fact 
very close to the p-mass obtained for the same action. 

Lattice calculations, although starting from first prin- 
ciples, are not free of approximations. The most obvious 
one is the use of the lattice itself, necessarily resulting 
in discretization errors. These errors can be reduced by 
the use of improved lattice QCD actions and the con- 
comitant improved observables. In the work reported 
here, we extend the previous work by working in 0(a) 
improved lattice QCD, which guarantees that errors in 
the matrix elements we extract are only of order 0(a 2 ). 
In order to emphasize the importance of consistently us- 
ing both improved action and observables, we discuss the 
form of the vector current operator at some length and 
give numerical examples for results one obtains with the 
current operators used in other work. 

In addition to the step from a discrete lattice to the 
physical continuum, one also has to extrapolate the lat- 
tice results in the pion mass. As is known, lattice calcula- 
tions yield results for pions much heavier than the ph ysi- 
cal pion. The previous form factor calculations in [l(llll| 
were for pions on the order of 1 GeV. Another improve- 
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ment step which we undertake in this paper is to extend 
our calculations of the electric form factor down to pion 
masses of 360 MeV. For the mean square charge radius of 
the pion which we extract from the form factor, we then 
study the extrapolation to lower masses. 

It is instructive to compare the mean square radii ob- 
tained from the Bethe-Salpeter amplitudes and from the 
form factor. We find considerable differences which, as 
suggested earlier @ can be ascribed to the effect of the 
gluon motion on the position of the center of mass. 

The paper is organized as follows. In Ch. 2, we first 
describe the general features of our approach and the de- 
tails of our lattice calculations. Results for the 2- and 
3-point function are described in Ch. 3 and 4, respec- 
tively. Our findings for the form factor and the mean 
square radius of the pion are elaborated on in Ch. 5. A 
summary of our work and conclusions are contained in 
Ch. 6. 



II. THE GENERAL METHOD 

A. The observables 

To obtain the pion form factor from the lattice, one 
has to calculate two observables, the 2- and the 3-point 
Green's function for an interacting quark- antiquark pair 
with the quantum numbers of the pion. Improved tech- 
niques are used in order to remove 0(a) discretization 
effects. We do this non-perturbatively for both the ac- 
tion and the electromagnetic vector current in the 3-point 
function. Fits to both observables were used to extract 
the desired information, such as the form factor and the 
pion mass. 

The 2-point function, shown schematically in Fig. ^ 
projected to momentum p is given by 




R 



FIG. 1: 2-point function. 
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The operator (ft creates a quark- antiquark pair with the 
quantum numbers of the pion at the source at (0,0), 
while (j>(x) annihilates it at the sink. Since we will con- 
sider a 7r + , it is given by 



4>\x) = ip u {x) 75 4> d {x) 



(2) 



Below, all flavor, spin and SU(3) color indices will be 
dropped. On the sink side, we use an extended oper- 
ator with an inter-quark distance R. This suppresses 




FIG. 2: 3-point function. 

the contribution of excited states to the 2-point function 
and facilitates the extraction of the pion mass. In or- 
der to keep the calculation gauge invariant, quark and 
anti-quark at the sink are connected by gauge links. To 
further enhance the contribution from the pion, the links 
in the extended pion operator have been fuzzed to better 
simulate the tube-like nature of the gluon cloud. 

The 3-point function, shown in Fig. [2 again concerns a 
quark-antiquark pair with pion quantum numbers, prop- 
agating from Xi to Xf ; disconnected diagrams do not con- 
tribute P, lllj | . At an intermediate time t a photon is 
coupled to either one of the charged quarks. This observ- 
able is obtained in momentum space by calculating 

X f ,x 

x £ -i Pf(x f -x) -i p,-x ^ ^ 

The parameter R for the pion operator at the sink 
is now fixed to the value giving the best overlap in the 
2-point function. As we will further discuss below, the 
choice of the current to which the quarks couple is im- 
portant. The continuum or local current, 



3% = iP(xh^(x) , 



(4) 



is not conserved on the lattice and needs renormalization 
by a factor Zy, yielding the renormalized local current 
jf} L . Using the Noether procedure, one can also con- 
struct a current which is conserved on the lattice 1131 . 



-$(x + (i)(l + 'Y li )Ul(x)i/>(x) 



(5) 



This conserved current still requires 0{a) discretization 
corrections for matrix elements away from the forward 
direction. 

Using Symanzik's improvement program, one can iden- 
tify appropriate operators, which, when used together 
with the improved action, result in matrix elements that 
are free of all 0(a) discretization errors. For the vector 
current considered here, the resulting improved current 

is Guam 



with 



il = z v {j%{x) + ac v d v T^} , 



T^u = ip(x) i cr^u ip{x) 
Zy = Z v (l + ab v m q 



(6) 



(7) 
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It is conserved to 0(a 2 ) and differs from the renormalized 
local current by a total divergence, which vanishes for 
forward matrix elements. The bare quark mass is defined 
as 



(8) 



where k c is the kappa value in the chiral limit and 
a the lattice spacing. The constants in j are non- 
perturbatively determined from lattice simulations and 
as such completely remove the 0(a) effects. 

From current conservation it can be shown that the 
general Lorentz structure of the matrix element for the 
electromagnetic current of an on-shell pion is 



{^{pi)\jA^(Pi)) c 



(p f + Pz ) fl F(Q 2 ), (9) 



where F(Q 2 ) with Q 2 = — (p/ — Pi) 2 is the form factor 
we are interested in. Connection with the continuum 
description is made by a proper normalization, 



{^iPf)\3nW(pi))l at 



(7r(P/)|jpk(P»)) f 



(10) 



where Ef and Ei are the final and initial energies, re- 
spectively. 

In our calculations of the 3-point function we project 
on initial and final three momenta with the same length, 

\Pi\ = Ip/I , (ii) 

which implies that there is no energy transfer to the pion, 

Ei = E } . (12) 

The 4- momentum transfer to the pion, Q 2 = (pf — Pi) 2 , 
is then varied by changing the angle between the two 
momenta. Since we will use the [i = 4 component of 
the current, this choice has, among others, the numerical 
advantage that we have 



Ef + E, 
2yJ EfEi 



1 



(13) 



when extracting the form factor F(Q 2 ). 

The Vector Meson Dominance (VMD) model has been 
quite successful in describing both experimental as well as 
early lattice data. This model is inspired by effective field 
theory and is schematically depicted in Fig.|Hl Assuming 
universality, i.e. g P irir — g P , the form factor is given by 
the simple monopole form 



' mi. 



(14) 



The lattice simulation 



We performed calculations in the quenched approxima- 
tion on a 24 3 x 32 lattice. A set of 100 gluon configura- 
tions at a coupling of (3 = 6.0 was generated. Thermaliza- 
tion was reached in 2500 sweeps, whereafter we obtained 
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FIG. 3: The VMD model. 



configurations at intervals of 500 sweeps. One sweep con- 
sists of a pseudo-heatbath step with FHKP updating in 
the SU(2) subgroups, followed by four over-relaxation 
steps. 

We used the improved Shcikholcslami-Wohlcrt ac- 
tion [Tfil | with the non-perturbatively determined p"?l | 
value of csw — 1-769. With this action, we computed 
propagators for five values of the hopping parameter cor- 
responding to pion masses 1 of 970, 780, 670, 540 and 
360 MeV, see Table 1. We imposed periodic bound- 
ary conditions except in the time direction where for 
the fermions anti-periodic boundary conditions were im- 
plemented. The values of the constants cy,bv and Z v 
needed to eliminate the 0{a) effects and to renormalize 
the current were taken from Bhattacharya et al. |l9j |. 



III. THE 2-POINT FUNCTION 

To extract physical information from the numerical 
data for the 2-point function, Eq. we use the following 
parametrization 



e p 2 



n=0 



x .•osi l {/-.;;::4 r H-i 



where N T = 32 is the extension of the lattice in the time 
direction. We include the contribution of the ground 
state (n = 0) with energy E® and of a first excited state 
(n = 1) with energy E^. As discussed in connection with 
Eq. ^ the parameter R indicates the quark-antiquark 
distance at the sink, which will be chosen to enhance the 
contribution from the pion ground state. For this we use 
the approach as originally proposed in [2(|. The fuzzed 
gluon links at the pion sink are created with a link/staple 
mixing of 2 and a fuzzing level of 4. The denote the 
matrix elements, 



Z£(p) = |<n|fo|n,p) 



(16) 



1 For definiteness we have taken a 
scale. 



: 0.105 fm from to set the 



4 



which also will yield the 'Bethe-Salpeter' amplitudes 
from which information about the structure of the pion 
can be extracted. 

The data for the 2-point function which correspond to 
the same absolute value of the spatial momentum are 
averaged per configuration. The different configurations 
are then combined in jackknife averages with a block- 
size ranging from 1 to 5. No significant changes in the 
errors of the averages were observed for increasing block- 
size, indicating that there are no significant correlations 
in our ensemble. 



A. The fuzzing distance R 

To determine the optimal value for the inter-quark dis- 
tance R we used the jackknife averages to calculate the 
effective energy of the pion, 

E eff (t, \p\) = ln«G 2| fl(*,p)) / < G 2,*(t + (17) 

We varied R from to 10, plotted the effective energy for 
these different fuzzing levels and looked which one stabi- 
lizes first. An example for the dependence on R is shown 
in Fig. 4(a) for the effective mass, M e 1 1 (t) = E e ff(t,0), 
and for the effective energy in Fig. |4(b)] with \p\ 2 — 0.48 
GeV 2 , corresponding to the momentum of the pion in 
our form factor extraction below. The optimal R-value is 
somewhat dependent on the pion's momentum and mass. 
After several such tests, we chose R — 3 as the common 
extension parameter for all calculations. Enhancing the 
ground state contribution is particularly important for 
the 3-point function, where the distance between source 
and sink is typically small. 
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B. Pion masses and energies 



FIG. 4: Effective masses (a) and effective energies (b) for 
various R. In (a) the pion mass is m w = 970 MeV, in (b) 
= 540 MeV and p 2 = 0.48 GeV 2 . 



Having chosen the fuzzing distance i?, we extract the 
pion masses and energies by fitting the jackknife averages 
to Eq. E| for both a single state and two states. The fit 
range, the number of included time slices centered around 
the midpoint of our time-grid at t = 16, is reduced until 
the minimum % 2 is found and consistency between both 
fits can be checked. In extracting the masses from the 
p = averages, we found that for a single-state fit a fit 
range of about 15-17 time slices gives the lowest \ 2 ■ I n 
case of a two-state fit, fitting the complete i-range (31 
time slices) yields the smallest statistical error; x 2 /dof 
is about the same for different fit ranges. We found con- 
sistency between the two fits. The resulting masses are 
given m TableUand plotted in Fig.0 as a function of the 
inverse of the hopping parameter k. They agree very 
well with the results obtained in e.g. [2lH2^ | who use the 
same action as we do. Also shown are extrapolations to 
the chiral limit, based on the following fit functions, 

ml = Cl ( - - —) (18) 



TABLE I: Masses and { r2 } BS f° r different k- values. 

k m q m n m p (r 2 ) flS 

0.13230 154 MeV 970(4) MeV 1188(6) MeV 0.1414(2) fm 2 

0.13330 101 MeV 780(4) MeV 1053(8) MeV 0.1480(2) fm 2 

0.13380 75 MeV 670(4) MeV 989(9) MeV 0.1508(2) fm 2 

0.13430 45 MeV 540(6) MeV 0.1526(2) fm 2 

0.13480 23 MeV 360(9) MeV 0.1528(4) fm 2 



ml = c 2 and (19) 

\K K C J 

~ = ^- + hml + b 2 ml , (20) 

resulting in a slightly different value for k c . Quenched 
chiral perturbation theory predicts the second form with 
5 small and positive • As in [2l| , we obtain a negative 
value for 5. Eq. [201 is a phenomenological fit [2l]]. There 
is not a significant difference in the fit quality, so we 
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FIG. 5: Pion masses as function of n. Lines: extrapolations 
as indicated. 
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The data deviate from both predictions at higher mo- 
menta. However, the figure clearly demonstrates that for 
the momentum relevant for our form factor extraction, 
p 2 = 0.48 GeV 2 , we are dealing essentially with contin- 
uum kinematics. 

C. Results for the p-mass 

In the discussion of the form factor, we will often refer 
to the vector meson dominance model. For completeness 
and as a further test for our methods, we have also ex- 
tracted the mass of the lowest vector meson, the p-meson. 
Proceeding analogously as for the pion, Eq.^ we consider 
a 2-point function with source and sink operators of the 
form 

Vi=$~firl>, (23) 

which project onto the polarization state i of a vector 
meson. With the same boundary conditions as for the 
pion, the 2-point function for three-momentum p = 
was then fitted to a cosh- form as in Eq. ^| We averaged 
the polarization states i — 1,2, 3. The results for the 
three highest parameters n are shown in Table [I] They 
agree with the values obtained in [2lT |22| . For the re- 
maining two k values our statistics were not sufficient. 
When needed later for the comparison with the vector 
meson dominance model for the form factor, we will for 
simplicity use the values from [2j. 

D. Bethe-Salpeter amplitude and ( r2 ) ss 



FIG. 6: Energy- momentum relation for k — 0.13230 and 
K = 0.13430; solid line: continuum relation, dashed: lattice 
dispersion relation. 



cannot prefer one fit over the other. Therefore, we simply 
average the different values, yielding n c — 0.13524(4), 
which agrees quite well with values obtained from the 
literature, k c = 0.13531(1) fH) and k c = 0.13525 pf . 
The bare quark masses we obtain with our k c are given in 
Tabled They will be used in the (small) mass dependent 
correction of the improved current. 

Proceeding analogously, we have also extracted the 
pion energies E® for several non-vanishing three- 
momenta, again using single and two-state fits in com- 
bination. The results for the ground state energy are 
shown in Fig. together with the prediction from the 
continuum dispersion relation, 



E cont = y/^l+p 2 , (21) 
and from the generally favored lattice dispersion relation, 



• i 2 & ■ 1,2 M \ - . 2 Pi 

sinn — = sinn \- > sin — 

2 2-^2 



(22) 



To obtain the 'Bethe-Salpeter (BS) amplitudes' or 
'wave functions' &(R), we use the Z factors extracted 
from a fit of the 2-point function (Eq. I15|) , since for a 
pion at rest, 



$(R) = JZ%(0)/Z°(0). 



(24) 



We simultaneously fitted the jackknife averages for R 
ranging from to 10. Again the results from both 
the single- and a two-state parametrization were used 
to check the consistency. 

The same results for the wave function, but with 
smaller errors, were obtained from a fit to the plateau 
of the ratio 



Hr) 



G a ,ji(t,0) 



(25) 



G 2 ,o(*,0) 

An example for a BS-wave function is shown in Fig.0for 
the heaviest pion; only the result of the second method 
is displayed. The wave function can be used to obtain 
an estimate of the mean square charge radius of the pion 
according to 0, 



1 Jd 3 ff 2 $ 2 (|r1) 
w JdV $ 2 (|r~j) 



(26) 
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FIG. 7: BS wave function for n = 0.13230. 

The factor w is included to reflect the uncertainty in the 
resulting (r 2 ). If one assumes that quark and antiquark 
are always located on opposite sides of the meson center 
of mass at distance r/2, one has w — 4; assuming that 
the quarks move uncorrelated around the center-of-mass, 
w = 2 should be used. The BS-radii as a function of the 
pion mass are given in Table [I] for w = 2. The values are 
much lower than the physical value of 0.439(8) fm 2 [24| . 
Moreover, the mean square radius is seen to be almost 
independent of the pion mass in the investigated range. 
We will comment on this in more detail in Section Ivl 



leaves the source in the forward and backward direction 
in time and guarantees that F(Q 2 ) — 1 at Q 2 = 0. This 
relation holds for each background gauge field configu- 
ration separately and thus also for configuration aver- 
ages. We have verified that our results for the conserved 
current satisfy this relation to an accuracy better than 

o(io- 5 ). 

In order to obtain F(0) for the renormalized local and 
the improved current, we again exploited Eq. [57| The 
lhs of this equation was averaged over pairs of values t 
and t' symmetric around tf and normalized by the 2- 
point function. Utilizing the Zy factor from [19j gives 
F J (0) = F RL (0) = 1 within a jackknifc error of 1 %. 
Alternatively, we could have applied this method to in- 
dependently extract Zy as done e.g. in _25j. However, 
for consistency we used the entire set of improvement 
parameters from |l9j| . 

While the above method allows one to determine F(Q), 
we now describe how we extract the form factor for 
Q 2 > 0. As in the 2-point function we allow two states to 
contribute and parametrize the 3-point function, Eq. |3 
as 

l l 

G 3 Atf,t;p f ,Pi) = E E \Jz%(Pf)Zo(Pi) 

m— n— 

x (m,p f \U0)\n, Pi ) e- E *f * , (28) 



IV. THE 3-POINT FUNCTION 

For the 3-point function, already introduced in Ch. 2, 
we now consider a pseudo-scalar source at t — 0, a sink at 
tf = 11 and a coupling to the photon at t with < t < tf. 
Barad et al. 0] pointed out that current conservation 
provides an important numerical test which relates the 2- 
and 3-point functions. Translated into momentum space 
and for our periodic boundary conditions, this relation 
reads 

G 3 (t f ,t;p,p) - G 3 (t f ,t';p,p) = G 2 {t f ,p) , (27) 

where, in the second term on the left hand side, tf < 
t' < N T . It amounts to considering the total charge that 



where (m, n) ^= (1,1). Contributions from, for example, 
the production of pion pairs, as well as 'wrap around 
effects' due to the propagation of states beyond tf are 
exponentially suppressed (< C(e -5 )); similarly, an elas- 
tic contribution from the excited state was estimated to 
be of the order of 1% or less. All these effects are not 
reflected in the chosen paramctrization. The inelastic 
transitions <-» 1 are included to better describe the 
data. However, it should be understood that the state 
labelled 1 parametrizes contributions from all possible 
excited states. Therefore we do not interpret our param- 
eters as the energy or the transition form factors corre- 
sponding to a single genuine excited state. 

Since for a given momentum the pion is the state with 
the lowest energy, one gets with Eas. l^landllUI 



G 3 , M (i/,t;p/,p l 



F(Q 2 ) ^±E^,/zUPf)Z°(p t )e-^ t 



+ V Z MP/)^o°(P,) (l,Pf\U0)\0,Pi) e 



(i-o) 



(29) 



In the simulations, we took \pf\ = \pi\ = v2|Pmin.|, where 



2ir 

\Pmin\ = -Tf— ( 30 ) 
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is the minimal momentum for a lattice with N a lattice 
points in the spatial extension. In our case, it amounts 
to 



Pi 



p) = 0.48 GeV 2 



(31) 



For our analysis we use the fourth component of the cur- 
rent, (i = 4. With our choice of momenta the kinematical 
factors in the fit function, Eq. 1291 therefore simplify con- 
siderably, see Eq. El note also that the i-dependence 
of the first term in Eq. EH vanishes. As a result, the 
form factor is more easily extracted without restricting 
the simulation too much. Different momentum transfers 
are obtained by varying the relative orientation of Pi and 
Pf- 
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A. Extraction of parameters 

We begin by averaging the 3-point correlation func- 
tions which have the same four-momentum transfer 
squared and then again combine the configurations in 
jackknife averages. Typical ja ckknif e averages of the 3- 
point function are shown in Fig. 8(a) for the one but high- 
est pion mass for various momentum transfers an d for d if- 
ferent masses at fixed momentum transfer in Fig. 8(b) If 



only the pion ground state would contribute, there would 
be no i-dependence in this quantity. The data however 
clearly indicate the admixture of an excited state. This 
is especially seen at high Q 2 and for low pion masses, 
where there is no time-slice where it is safe to assume 
that only the pion ground state is present. 

We therefore chose to proceed by simultaneously fit- 
ting the parameters in the 2-point and 3-point function, 
Eq. El an d Eq. |2H1 respectively. We hereby exploit the 
fact that certain parameters appear in both Green's func- 
tions. In the case of the 2-point function, we fitted the 
energies, E° and E 1 and the Z factors over the complete 
time interval, 1 < t < N T — 1. In the 3-point function, 
we fit the same energies and Z-factors, and in addition 
the form factor F(Q 2 ) and the transition matrix elements 
over the interval ti < t < tt. We assume that the ener- 
gies and Z factors only depend on the magnitude of the 
three-momenta and use the fact that we chose | p/ 1 = | Pi | . 
The fits are done for each value of Q 2 separately. The 
values for x 2 /dof lie between 0.15 and 0.40, depending 
on mass and momentum transfer. The energies and Z 
factors we obtain from our fits at different Q 2 agree to 
high accuracy because they are largely determined by the 
2-point function. 

To compare with earlier work |lOtlll| we also extracted 
an estimate for the form factor from the ratio of 3- and 2- 
point functions. However, the assumption of just a single 
state contributing is at the basis of this method. Corre- 
spondingly, we found differences ranging from 5 — 10% 
for F(Q 2 ) between the ratio method and our combined 
fit method, where inclusion of an excited state clearly im- 
proved the fit quality. The size of the difference depends 
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FIG. 8: Improved 3-point function for (a) different Q 2 
m n = 780 MeV and (b) different rrw at Q 2 = 0.97 GeV 2 . 
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on the pion mass and the momentum transfer, which in- 
fluences the flatness of the 3-point function in the middle 
between source and sink. All our results in the next sec- 
tion are therefore based on the fit method. 



V. RESULTS AND DISCUSSION 

We will now discuss the form factors obtained from the 
procedure described in the previous Chapter. 

A. Dependence on the current 

We first compare the form factors we obtain with the 
local current, Eq. 0] the conserved current, Eq. [5J and 
the improved current, Eq. Only the improved current 
ensures that there are no corrections to 0(a) with our 
action. 

The form factors for Q 2 > 0, obtained through the 
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FIG. 9: Form factors extracted from different currents as a 
function of Q 2 for 2 pion masses: m n — 970 MeV (top) and 
= 360 MeV (bottom). Solid curve: VMD model predic- 
tion with my = m p taken from 2l| . Data shifted horizontally 
for clarity. 
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FIG. 10: Improved form factor as a function of Q 2 for differ- 
ent m-K. Lines: fits to VMD-form, Eg. 1141 

to the result for the renormalized local current, also for 
the light pion and out to our largest momentum transfer. 
That these two form factors are so close is due to the fact 
that the contribution of the tensor term in the improved 
current is small. Since the matrix element of the tensor 
operator can become comparable in size (up to 70 %) to 
the local current operator, this smallness is due to the 
fact that the coefficient, cy = —0.107, determined in 0] 
is rather small; the preliminary value obtained by the 
ALPHA collaboration is much larger, cy = —0.32. 
Since the improved current is a linear combination of the 
local and the tensor term, Eq. the change in the im- 
proved current due to a change in cy is straightforward. 
The difference between both values of cy is of order a, 
resulting in improved currents which are different only at 
order a 2 . This shows that 0(a 2 ) effects can still become 
as large as 10 % at higher Q 2 values and low masses. 

B. The form factor and vector meson dominance 



simultaneous fit procedure, are shown in Figs. 9(a) 
and 9(b) for two different masses. As can be expected 



from Fig. 8(a) the same fit procedure yields values for 
F(0) with an error comparable to the results for low Q . 
However, we use the method discussed in connection with 
Eq. [23 to extract F(Q) — 1 to higher accuracy. 

We first of all observe that, as expected, the results for 
our heaviest pion, = 970 MeV, are much more ac- 
curate than for the lightest pion, m, = 360 MeV. How- 
ever, even with the larger error bars, a rather smooth 
Q 2 -dependence is seen also in the latter case. We fur- 
ther observe that the differences between conserved and 
improved current grow with momentum transfer and de- 
creasing mass, resulting, in particular for the light pion 
mass, in a substantial correction. The differences be- 
tween the improved and the renormalized local current 
are due to the tensor term. Although this contribution in- 
creases with Q 2 the improved form factor stays very close 



As was already observed in [Tl|], the lattice results 
for the form factor show a behavior expected from vec- 
tor meson dominance. In Figs. |9(aJ| and 9(b) we show 
the prediction for the form factor if we use the simple 



monopole form, Eq. 1141 with my 



the lattice p- 



mass at the same k- value |2jj. At large pion mass the 
VMD-prediction describes the form factor based on the 
conserved current rather well but lies substantially above 
the results from the improved current. However, at lower 
pion masses the model prediction shifts toward the im- 
proved form factor results. To investigate this point in 
more detail, we fitted the improved form factors, using 
the vector meson mass my as a fit parameter, omitting 
the point at the highest Q 2 value. The parametrization 
works well and results are shown in Fig. ^] Table [H] 
compares the fitted my values to the p-mass extracted 
from 2-point functions j2l|. The two values can be seen 
to come closer together as the pion becomes lighter, sug- 
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TABLE II: m p from the 2-point function [2l| . fitted my and 
(r 1 ^ for different values. 



970(4) MeV 
780(4) MeV 
670(4) MeV 
540(6) MeV 
360(9) MeV 



1169(3) MeV 
1032(4) MeV 
966(6) MeV 
901(6) MeV 
841(24) MeV 



my 

1086(26) MeV 
968(26) MeV 
931(26) MeV 
882(36) MeV 
833(75) MeV 



0.197(9) fm 2 
0.249(13) fm 2 
0.269(15) fm 2 
0.299(24) fm 2 
0.34(6) fm 2 



gesting a better agreement of the improved results with 
the simple vector meson dominance model for lower pion 
masses. However, as we will see later, in the physical 



limit, using my 
data accurately. 



fails to describe the experimental 



C. Determination of the charge radius 

It is well known that the slope of the form factor is 
related to the mean-square charge radius of the pion, 



dF(Q 2 ) 



dQ 2 



Q 2 =0 



(32) 



In contrast to the charge-radius extracted from the 
Bethe-Salpeter amplitude, this determination of (r 2 ) is 
not based on any specific assumptions about the quark 
motion inside the pion. As the vector meson dominance 
model works very well for low Q 2 , we show for simplicity 
the mean square charge radius of the pion obtained from 
the monopole fit, which yields 



(r 2 ) 



G 



(33) 



In the following we only work with the improved cur- 
rent. By looking at the values in Table ITT1 and in Fig. II II 
we observe that the (r 2 ) extracted from the form factor 
shows a considerable mass dependence. This is in con- 
trast to the BS results, which are also shown. Moreover, 
these results, which we obtained in Section l!!! Dl are con- 
siderably lower than the value we extract from the form 
factor. As already discussed by Gupta et al. |6j, this can 
be due to how one treats the center of mass of the two 
quarks. However, as these authors also point out, the 
value extracted through the form factor contains contri- 
butions that are not included when calculating (r 2 ^ from 
the Bethe-Salpeter amplitude. 

The (r 2 ) values obtained from the form factor can be 
seen to get closer to the physical value of (r 2 ) = 0.439(8) 
fm 2 as the pion mass decreases. This led us to try sev- 
eral extrapolations to the physical limit which will be 
described in the next section. 

In addition to the two methods discussed above, there 
is another method to obtain the charge radius of the pion 
from lattice QCD. This method is based on calculating 
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FIG. 11: Pion charge radii extracted from the form factor 
and extrapolations in m 2 . The BS-radii are included for com- 
parison. Experimental result from |24| 



density-density correlations or 4-point functions for the 
pion [ij, |8j . It has recently been used by Alexandrou et 
al. Q for densities at equal times. However, there are 
difficulties in the extraction of (r 2 ) from density-density 
correlations as discussed in detail by Burkardt et al. 
and Wilcox [27| . 



D. Extrapolation in m n 

To obtain more physical results, one can try to extrap- 
olate in the pion mass. We take (r 2 ) as the quantity 
to extrapolate, since it is known experimentally and its 
extrapolation has been discussed in the literature. We 
consider three different types of extrapolations. From 
chiral perturbation theory (x_PT), one knows the 1-loop 
result lH j 



( 



2 \ one— loop 
r 'XPT 



ci + c 2 In m 2 . 



(34) 



In our fit, we will treat ci and c 2 as free parameters. In 
quenched yPT, the radius is constant at this order of 
expansion 29] . It is however expected that this situation 
will change at the two-loop level, which will introduce 
terms like 



1 



di In m„ 



d 3 ml , 



(35) 



including a term linear in m 2 which, for our pion masses, 
can be expected to yield the dominant contribution |30| . 
We therefore only tried a form containing a constant plus 
a term linear in m 2 

We have observed that our form factor data can be 
well described by a monopole form as suggested by simple 
vector meson dominance. Therefore, we use this model 
to obtain an additional extrapolation. Since it can be 
seen that my, just like m p , scales approximately linear 
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FIG. 12: Improved form factor for different m^, compared to experiment [3ll l32| : broken lines as in Fig. 1101 



with m 2 , one arrives at, 



6 



VMD 



(bi + b 2 ml 



1 \2 



The three extrapolations are plotted in Fig. ^] The 
extrapolated value is seen to depend strongly on the 
method chosen. The VMD ansatz describes the data 
best. If we use this ansatz to extrapolate to the phys- 
ical pion mass we find (r 2 ) = 0.37(2) fm 2 , which lies be- 
low the experimental value of 0.439(8) fm 2 . This could 
clearly be due to the extrapolation chosen, but also due 
to assumptions and approximations, such as quenching, 
inherent in our approach. 



between our lattice calculations and experiment is not 
clear. The solid line in the figure shows the monopolc 



(36) form using rn\ 



= 6/(r 2 ) 

' \ / exp 

;r 2 ) 



with the experimentally 
measured charge radius (r 2 ) . This describes the ex- 

° \ / exp 

perimental data quite well also away from Q 2 — 0. How- 
ever, the corresponding vector meson mass at 730 MeV 
is significantly lower than the p-mass, emphasizing that 
the VMD-inspired monopole description provides a suc- 
cessful parametrization of the form factor data, but does 
not hold in detail. 



VI. SUMMARY AND CONCLUSIONS 



E. Comparison with experiment 

In Fig. ^] we show our results together with the avail- 
able measurements |33| • For clarity we only show our 
results for three K-values. As can be clearly seen, all our 
calculated form factors lie above the measured values. 
Nevertheless, a continuous trend toward the experimen- 
tal values can be observed and we come quite close to 
them. Whether only a straightforward further lowering 
of the pion mass will resolve the remaining discrepancy 



We have presented calculations for the pion form fac- 
tor that improved and extended previous work in sev- 
eral respects. We have pushed the form factor calcula- 
tions for a large range of Q 2 towards much lower pion 
masses than before. In doing this, we have worked in a 
framework that guarantees the absence of 0(a) correc- 
tions. This meant a consistent use of an improved action 
with the concomitant improved conserved current. It was 
shown that use of this improved current leads to signifi- 
cant changes over results based on the conserved Noether 
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current for this action, which still contains 0(a) correc- 
tions at finite Q 2 . We chose kinematics where the initial 
and final pion momentum had the same absolute value, 
which leads to practical simplifications when extracting 
the form factor. For the momenta we use, we confirmed 
that energy and momenta are sufficiently close to satis- 
fying a continuum dispersion relation. 

Our results for the form factor were seen to smoothly 
vary with pion mass. The lower Q 2 results could all be 
described quite well by a simple monopole form factor. 
The fitted range parameter rriy 1 was, for each K-value, 
found to be close to the corresponding lattice p-mass. 
The agreement between the two values got closer for de- 
creasing pion mass, indicating better agreement with the 
vector meson dominance model. 

The form factor can be used to extract the mean square 
charge radius of the pion. The values we obtained show 
that the estimates for (r 2 ) based on the Bethe-Salpeter 
amplitude are qualitative as well as quantitative not very 
reliable. Disagreement of up to a factor two was found 
with the form factor based values, which showed also 
a quite pronounced mass-dependence in contrast to the 
BS results. Extrapolations of our charge radii towards 
the physical pion mass were shown to lead to no unique 
prediction. The best description of the results at our 
pion masses was provided by a vector meson dominance 
model. When extrapolated to the physical pion mass, it 
yields a value for (r 2 ) about 15% below the experimental 
value. For an extrapolation inspired by (quenched) chiral 
perturbation theory our pion masses are too high to be 
sufficiently sensitive to the predicted lnmj terms. 

When compared to the experimental form factor, it 
could be seen that our results consistently approach the 
data from above over the entire range of Q 2 we consider. 
While gauge invariance fixes all form factors at Q 2 = 
to F(0) — 1, we see that the calculated form factor at 
Q 2 > comes close to the experimentally determined 
shape, and to a monopole parametrisation. This is a 
nice confirmation that lattice QCD indeed describes a 
non-perturbative feature such as a pion form factor quite 
realistically and in detail. However, a straightforward 
extension of our approach to even lower pion masses or 
higher Q 2 is not necessarily the way to proceed to close 
the last gap to the experiment. Improvements of our 
approach and other lattice methods will become neces- 
sary. Corrections of order 0(a 2 ), for example, will be- 
come increasingly important and one also has to under- 



stand the role of dynamical quarks, which are neglected 
in the quenched approximation. 

As is well known, Wilson fermions have the major dis- 
advantage that chiral symmetry is broken, already at 
0(a). This was one of the reasons improvement was in- 
vented and why we chose a framework where action and 
operators where consistently improved and only correc- 
tions to order 0(a 2 ) show up. Another method, among 
others, is the introduction of a fifth dimension and use of 
so-called domain wall fermions. Chiral symmetry is then 
not tied to taking the continuum limit. The price one 
pays is a substantial increase in the computer time. The 
RBC collaboration has chosen this approach and first re- 
sults can be found in |33| . In this paper pion masses down 
to 390 MeV are used, albeit on a coarser lattice. Their 
results obtained so far at two low Q 2 points, based on 
the renormalised local current, seem to agree reasonably 
well with our values. Differences in the implementation 
of chiral symmetry show up at 0(a 2 ). 

An open question is of course the significance of the 
quenched approximation. Alexandrou et al. have cal- 
culated density-density correlations for the 7r, p, N and 
A in quenched as well as unquenched lattice QCD. In 
contrast to the p and A, only rather small effects are 
seen for the tt for values of around 600 MeV. The 
study of effects due to dynamical quarks, clearly more 
important at lower pion masses and high Q 2 , is an area 
where further work is necessary. 
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